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Abstract: 

We investigate the dynamics of an iV-level quantum system weakly coupled 
to a thermal reservoir. For any fixed temperature of the bath there exists 
a natural reference state: the equilibrium state of the system. Among all 
quantum operations on the system one distinguishes Davies maps, they pre- 
serve the equilibrium state, satisfy the detailed balance condition and belong 
to a semi- group. A complete characterization of the three dimensional set 
of qubit Davies maps is given. We analyze these maps and find their mini- 
mum output entropy. A characterization of Davies maps for qutrits is also 
provided. 

PACS: 02.10.Ud (Mathematical methods in physics, Linear algebra), 03.67.- 
a (Quantum mechanics, field theories, and special relativity, Quantum infor- 
mation), 03.65.Yz (Decoherence; open systems; quantum statistical meth- 
ods) 

1 Introduction 



The quantum dynamics of simple model systems is a subject of consider- 
able research interest in view of applications in quantum information pro- 
cessing [H |2]. Any discrete dynamics of an iV-level system, which can be 
realized in a physical experiment, can be described by a quantum operation, 
also called stochastic map. This is an affine transformation of the set of 
quantum states of size N that is completely positive, see e.g. OH]. 
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Within the set of quantum operations one distinguishes the subset of bis- 
tochastic maps, they additionally satisfy the unitality condition: the uniform 
state is preserved. On one hand the set of bistochastic maps is easier to study 
and these maps enjoy several particular properties, for instance, the minimal 
output entropy is known to be additive for one-qubit bistochastic maps [5] . 
On the other hand, maps that are not bistochastic do play a crucial role in 
the description of several physical problems. 

Even in the case N = 2 the set of one-qubit operations is 12 dimensional 
and it grows as iV 4 for large N. Thus it is legitimate to ask, which of these 
maps could be of any special importance. To suggest any concrete class of 
maps one needs to make some physical assumptions and to investigate their 
implications. 

The main aim of this work is to select a physically motivated class of quan- 
tum operations, which includes also non-unital maps, and to study their 
properties. In theoretical investigations one often deals with isolated quan- 
tum systems, which are entirely separated from the environment. Such an 
assumption is not realistic, so we are going to assume that the iV-level sys- 
tem under investigation is weakly interacting with a thermal reservoir, see 



Figure 1: An N- level system weakly coupled to a reservoir at temperature 
T, 

The coupling with the bath at temperature T determines the equilibrium 
state pp of the system where the subscript (3 = 1/kT denotes inverse temper- 
ature and where k is Boltzmann's constant. A discrete quantum operation 
$ acting on the system will be compatible with the environment if it pre- 
serves pp. We assume, moreover, that the map arises from a continuous 
time evolution starting at the identity and that it satisfies the quantum de- 
tailed balance condition [Bj. These assumptions define the Davies maps [7j 
investigated in this work. 

The paper is organized as follows. In Section II we review some basic prop- 
erties of quantum states and quantum maps, discuss the detailed balance 
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N -level system 




reservoir 
in temperature T 
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condition and formally introduce the notion of Davies maps. In Section III 
we provide a complete characterization of the one-qubit Davies maps, dis- 
cuss possible uses of these maps in quantum optics and study their minimum 
output entropy. The case of Davies maps for three level systems is treated 
in Section IV. 

2 Definitions and conventions 
2.1 Qudits and quantum maps 

The states of an TV-level quantum system are described by N x N density 
matrices p, these are Hermitian, positive semi-definite and have trace one 

p = p\ p > 0, and Tr p = 1. 

States model the statistical data obtained by measuring observables A, in 
case of Hermitian matrices of dimension N 

(f(A))=Tt{pf(A)). 

Furtheron A4 at will denote the algebra of N x N matrices with complex 
entries. 

A quantum operation or quantum channel 3> is an affine transformation of 
the density matrices 

p^p' = Hp) (2.1) 

which is completely positive (CP). Remark that a quantum operation can be 
uniquely extended to a trace preserving linear transformation of the N x N 
matrices. Complete positivity means that the trivial extension $ (g) id^ to 
a space of any dimension d preserves positivity. The state transformation 
rule (|2.1|) is the Schrodinger picture of the action of the channel. One can use 
the corresponding Heisenberg picture as well, the map will still be completely 
positive but trace preserving becomes unity preserving. 

It was shown by Kraus [8] that a linear map <3? is completely positive if and 
only if it is of the form 

*{p) = Y t K i pK\. (2.2) 

i 

If $ is trace preserving, i.e. if $ is a quantum operation, then the Kraus 
operators K{ satisfy the relation 

i 
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There is a useful condition that guarantees that a linear transformation <5 
of the N x N matrices is completely positive. It suffices to check whether 
the map $ <g> id^ transforms a maximally entangled state 



1 N 



(2.3) 



into a positive matrix [9j [TU] . For a quantum operation $ 

1 - 

($®id JV )(|* + )(* + |) = ^ E $ (N)0'l)<»N)0'l>o (2.4) 

is a density matrix called the Jamiolkowski state or Choi matrix of <£. 

A quantum map $ is often represented in superoperator form $5. In this 
parametrization a density matrix p is written as a vector p by ordering its 
elements in a column. Equation (|2.ip can now be written as the action of a 
matrix &g on p 



$ s p. 



(2.5) 



For a single qubit one can take the basis consisting in the identity and the 
three Pauli matrices a = {<7i, 02, 0-3} = {a x , a y , a z }. A qubit state can then 
uniquely be written in terms of a vector in the closed centred unit ball in 
three dimensional real space, called Bloch ball 



P 



i (id+x-a), x & R 3 , ||x|| < 1. 



The quantum operation $ in this basis will be called a quantum operation 
in the Bloch form &b- As $ is affine, it transforms the Bloch ball into an 
ellipsoid inside the ball. Up to two orthogonal rotations of the basis the 
Bloch form of a one-qubit operation can be written as 
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The vector k describes the translation of the centre of the ellipsoid with 
respect to the centre of the Bloch ball and is zero for a bistochastic map. 
The absolute values of the rji are the lengths of the axes of the ellipsoid. 
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2.2 Quantum Markov process 



The reduced evolution of a quantum system coupled to a general environ- 
ment is hard to describe. Some additional assumptions make the problem 
more tractable. For example Markovianity means that the system interacts 
so weakly with the environment that it preserves its identity and that the 
evolution of the system does not depend on a whole history but only on its 
present state. A Markovian evolution in Schrodinger picture is described by 
a differential equation and an initial condition 

^-p = Qp and p(0) = p (2.6) 

where Q is called the generator. It consists of a unitary part 5 and a dissi- 
pative part L 

G = i5 + C. (2.7) 

Here 

5:p^[p,H] (2.8) 

is a commutator with an Hermitian system Hamiltonian H and C is a Lind- 
blad operator 

C : p _> J2(K a pKt - \ {KlK a , p}) . (2.9) 

a 

The K a S M.n can be compared with the Kraus operators (|2.2j) . The 
splitting (12. 7p into a unitary and a dissipative part is not unique. The 
solution of equation (|2.6p is 

p{t) = e gt Po . 

The form of a Lindblad operator (|2.9p guarantees that for any time t > the 
map exp(C?i) is completely positive, as was shown by Gorrini, Kossakowski, 
Sudarshan [Tl] and Lindblad |12j . 

A linear transformation Q on the N x N matrices is a generator of a 
semi-group of CP maps if and only if it satisfies a criterion similar to the 
Choi-Jamiolkowski criterion (12. 4p for CP maps. For Q that kill the trace, 
TrQ(A) = for all A, one has to check the positivity of the Choi matrix 
of the generator on the subspace perpendicular to the maximally entangled 
vector (|2.3p . Such a Q is then decomposable in the form (|2.9p . This type of 
condition is called conditional complete positivity [15]. 
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2.3 Quantum detailed balance condition 



Our model is an iV-level quantum system weakly coupled to a thermal bath 
at temperature T. The quantum system is described by its effective Hamilto- 
nian H which is diagonal with eigenvalues {e{\. The coupling with the bath 
distinguishes the equilibrium state pp at inverse temperature (3 = 1/kT, 
where k is Boltzmann's constant. The density matrix pa is called a Gibbs 
state, it is also diagonal in the eigenbasis of the Hamiltonian 



The weak-coupling limit yields semi-groups in continuous time of completely 
positive maps which enjoy the quantum detailed balance property, for a de- 
tailed treatment see [13] . Such semi-groups are called Davies semi-groups 
and they are to be compared with classical Glauber dynamics. The gener- 
ator of a Davies semi-group splits into a unitary and a dissipative part as 
in (|2,7p in such a way that the unitary and the dissipative parts decouple. 
Moreover, the dissipative part causes jumps between the eigenstates of the 
system Hamiltonian with transition rates that obey the micro-reversibility 
condition: the particle flow in equilibrium from level % to level j is exactly 
balanced by the reverse flow. Furthermore, off-diagonal elements of a density 
matrix decay at rates that are suitably adapted to the diagonal transition 
rates. A more precise definition follows below. 

A Davies map is an element of a Davies semi-group at some particular time. 
Due to the detailed balance property of the generator, a Davies map will sat- 
isfy the quantum detailed balance property with respect to the equilibrium 
state pfs [6]. This property is stronger than simple invariance and is usually 
expressed in terms of a self-adjointness condition with respect to the scalar 
product defined by the equilibrium state. The converse is, however, not 
true: there are quantum detailed balance maps that don't sit on a Davies 
semi-group, for some general results on the characterizing maps that sit on 
semi-groups see |15| . 

A Davies map is defined as follows 



Pp=2 exp(-/?#) or {pp)^ = — e ^ <%. 
Here Z is a normalization factor called partition function 



Z = 




1. $ = exp(£/i) 
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2. g = iS + £ with 5 and £ as in (J23J and (g^D , 



3. 5 and £ commute, and 

4. £ is self-adjoint with respect to the Gibbs state defined by the unitary- 
dynamics of the system in the following sense, for all A, B £ Mn 

Tr £ (B pp) A = Ti £(p/3 A) B. (2.10) 

Let us consider the matrix representation of the dissipative part £. The 
eigenvectors of the derivation 5(-) = [H , •] are the elementary matrix units 
with corresponding eigenvalues ej — ej. The third property from the 
above list guarantees that £ has the same eigenspaces as 5. Therefore the 
diagonal elements of a density matrix are not mixed with off-diagonal ones 
during the evolution determined by £. Furthermore, if the energy differences 
in the Hamiltonian are not degenerated off-diagonal elements are not mixed 
among themselves. Therefore the superoperator matrix £$ is block diagonal 
in the decomposition 

M N = Span({|i)(i| | i = 1,2, . . . N}) ©(© C\i)(j\) . 

The fourth property guarantees that the restriction of £ to the diagonal ma- 
trices is a classical detailed balance generator with respect to the probability 
vector of the eigenvalues of the Gibbs state and that the eigenvalues of £ 
corresponding to non-diagonal matrix elements are real. These conclusions 
also hold true for a quantum operation $ = exp(^t) and therefore for a 
Davies map. The requirement that Q and hence $ should be completely 
positive will be imposed by checking the positivity of the Choi matrix of the 
generator on the subspace perpendicular to the maximally entanglement 
state. The unitary evolution generated by 5 will play no role in what follows 
and we therefore simply omit it. It follows also from the fourth condition 
that the equilibrium state pa is invariant. 

5. Often one imposes ergodicity which requires that pp is the unique 
invariant state: for any density matrix p 

lim exp(Qt)(p) = lim exp(£t)(p) = pa 

t— *oo t—>oo 

Consider the following scalar product on the observables 

(X,Y) p :=Tr{p- l XiY). 
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The detailed balance relation (|2.10p for the dissipative part of the generator 
can be rewritten as 

(X,C(Y)) = (C(X),Y) p . 

Hence C is Hermitian for this scalar product and the same is true for exp(£i). 
As a consequence both C and exp(Ct) have a real spectrum. 



3 Davies maps on a two-level system 
3.1 One-qubit maps 

We write the Hamiltonian of the system in its eigenbasis 

H = \ n 1 I With e l > e 2- 

\0 e 2 J 

This implies that the Gibbs state consisting of elements of a probability 
vector (p, 1 — p) is of the form 

n _ (p \ _ 1 (expi-Pei) \ 

Pp -\P l-p)-Z\ exp(-/3e 2 )J' (3 - U) 

where (3 = 1/kT and the partition function reads 

Z := exp(— (3e\) + exp(— /3e2). 

The general superoperator form of a qubit Davies map is 



(I -a ap/(l-p) \ 

c 

0c 

\ a 1- ap/(l-p)J 



(3.12) 



where a and c are parameters of the map, while (p, 1 — p) is the vector 
of diagonal entries of the stationary state (|3.1ip . At infinite temperature 
P = the map is bistochastic, i.e. the uniform state is stationary. We still 
need to impose the condition that $5 lies on a one-qubit Davies semi-group 
of completely positive maps. This provides a full characterization of the 
allowed values of the parameters a and c, 



0<a<l, a + p<l, and < c < - a/(l -p). (3.13) 



S 



Because of the relation = exp(£i) the parameters a and c change in time 
according to 

a = (1 — p)(l — exp(— At)) and c = exp(— Ti), 
where A and T are constants such that 

r>ii>o. 

This constraint on the decay rate of off-diagonal elements with respect to 
the rate at which equilibrium is attained has been known for a long time, see 
e.g. |13} I14j. The constants A, T and the stationary vector (|3.1ip determine 
a generator £ of a semi-group {exp(Ct) \ t > 0}. 




Figure 2: Panel a) Allowed range of parameters in the (c, a) plane for one- 
qubit Davies maps (|3.12p shown for three different temperatures: T = 0, 
T p £ (0,oo), and T — > oo. The corners of the plot denote the identity 
channel 1, coarse graining ^cg an d a decaying channel $o : P ~ > |0)(0|. 
Panel b) The dashed lines in the allowed region at temperature T p represent 
the continuous time evolution exp(£t) for two typical generators C\ and £2- 

The allowed values of the parameters a and c occupy the region shown 
in Fig. [2] which depends on temperature. For higher temperatures the al- 
lowed region shrinks. In the limiting case T — > 00 the equilibrium state is 
maximally mixed, p/3=o = \ 1. This case corresponds to unital and thus 
bistochastic maps. The set of all one-qubit bistochastic maps forms a reg- 
ular tetrahedron spanned by the three Pauli channels a%, 02, and 03 and 
by the identity channel oq. The set of bistochastic maps lying on a semi- 
group is a proper subset: two tetrahedra glued at their base and localized 
in the vicinity of the identity map. The bistochastic Davies maps form a 
two-dimensional sub-region, see Fig. [3l 
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Figure 3: a) The tetrahedron of one-qubit bistochastic maps contains the 
shaded region P of maps of the form exp(Ct). b) The bistochastic Davies 
maps belong to the cross section of the set P with the plane determined by 
o"o = 1, 03 = <j z and the completely depolarizing channel <I>cd : P l— * \ 1- Its 
boundary is given by (13.130 , 



In Bloch parametrization a qubit Davies map takes the form 
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1 - a/(l - 


p), 


and 


^3 = 



where r]i = 772 = c, 773 
map acts on the Bloch ball as shown in Fig. [H 



(3.14) 



a(2p - 1)/(1 -p). Such a 



3.2 Davies maps in quantum optics 

Quantum dynamics consisting of a unitary part and a dissipative part are 
often considered in optics, the unitary evolution is caused by a driving exter- 
nal field while the dissipation is a consequence of a coupling to an external 
reservoir. The simplest instance is that of a two level atom. Its states are 
determined by a three component Bloch vector x = (u,v,w). The compo- 
nent w, the difference between the diagonal elements of the density matrix, 
measures the inversion of populations and the off-diagonal elements u and 
v are associated with an atomic dipole operator, u in phase and v out of 
phase with respect to the external field. The dissipation due to the inter- 
action with the environment decreases the length of the Bloch vector and 
therefore drives a pure initial state towards a more mixed state. 
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w 




Figure 4: The Bloch ball is transformed by the Davies map (|3.14p into an 
ellipsoid centred at x = (0, 0, w eq := 2p — 1) with axes (771 , 771, 773). 

The total evolution which is a composition of a unitary and a dissipative 
part is described by the Bloch equation, this is an affine differential equation 
for the Bloch vector. Affinity is a consequence of the linearity of quantum 
mechanics: an initial mixture of pure states evolves to the same mixture of 
evolved pure states. A few empirical constants t u , t v , t w , A, fi, and w eq 
enter the equations for x = (u, v, w) 



r d 


1 


— u = 




dt 


T u 


d 


1 


— v = 




dt 


T v 


d 


1 


—w = 




I dt 


T<w 



- Av 
+ An + 

{w — Weq) — 



The constants have a physical interpretation: t w is the rate at which inver- 
sion of the populations tends to the equilibrium state, t u and t v are decay 
rates of the atomic dipole, A is the field-atom detuning, and O is the Rabi 
frequency of oscillation between the ground and excited state of the atom. 
The parameter w eq fixes the equilibrium inversion. 

A semi-group of completely positive maps yields analogous equations for the 
Bloch vector. The generator £e of a Davies map in Bloch parametrization 
is a 4 dimensional matrix acting on the vector (1, u, v, w) T where (u, v, w) is 
the Bloch vector. We use the notation t\ = t% = t u = t v and T3 = t w . The 
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relaxation times t\ and 73 and the equilibrium inversion u> eq determined by 
an energy difference and a temperature are physical parameters while 771, 
772, and K3 describing the ellipsoid obtained as the image of the pure qubit 
states is rather a geometrical parametrization 
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Wq/T3 








-i/W 



(3.15) 



The variable w eq is related to the equilibrium state pp 
w, 



P 



■'eq 





l-p 

2p — 1. The Bloch equations describing a Davies map are: 



by 



r d 
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— u = 




dt 


T\ 
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—v = 
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—w = 




I dt 
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(w 



w, 



eq; 



By changing the basis one obtains the generator matrix in superoperator 
form 

(-{i-p) In jVt 3 \ 

-1/ti 

-1/ti 

\ {l-p) /n - p /t 3 J 



It can be checked that the outer (l,4)-block of this matrix satisfies the 
classical detailed balance condition with respect to the invariant probability 
vector (p, 1 — p). 

The quantum Davies map exp(£t) reads in Bloch parametrization 



exp(C B t) 



/ 


1 








° \ 























e-*/n 





V(i 


-e-*/ T 3) Weq 








e -*M / 



(3.16) 



The diagonal entries of this matrix belong to the spectrum of exp(Ct) which 
is real. In superoperator form (|2.5p the same map is represented by the 
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(l- e -*/*s)p \ 



1 - (1 -e-*/ T3 )p/ 

The conditions for lying on a semi-group of completely positive maps are 

Tl = r 2 >0, r 3 >0, andri<2r 3 . (3.17) 

These conditions are known as necessary conditions for a master equation 
of Lindblad type [16J. In this article the author obtains also a map of the 
form (|3.15p and the relations (|3.17p . In [17\ the map (|3.16p was considered 
in connection with the squeezed vacuum channel. It was derived directly 
from a specific form of the Lindblad generator. 

3.3 Minimal output entropy 

There are several possibilities for quantifying the randomizing properties of 
a channel both classical and quantum. We focus here on the minimal von 
Neumann output entropy |20U21j but Renyi entropies [221123] are also consid- 
ered. A different approach, the map entropy [24\ I25j. uses the Jamiolkowski 
state and, at least for classical channels, the entropy of an associated Markov 
process has also been investigated [15] . 

The minimal output entropy is defined as 

S min ($) :=minS($(p)), (3.18) 
p 

where S is the von Neumann entropy 

S(p) := — Tr(p log p), p a density matrix. 

Because the entropy is concave we can always find a pure minimizer for (|3.18p . 
In Fig. U] an inverse image of a point on the ellipsoid nearest to the Bloch 
sphere is such a pure minimizer. As was shown before the image of the 
Bloch ball under a Davies map has rotational symmetry with respect to 
an axis of the Bloch ball. Therefore, in this case, there is whole family of 
minimizers. In order to calculate the minimal output entropy it suffices to 
consider density matrices with all entries positive. This generalizes to higher 
dimensional problems. 
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exp(£ s t) 
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In a number of cases the minimizer for the output entropy is known. The 
simplest situation is probably a classical N state stochastic map. This is 
given in terms of a stochastic matrix: a N x N matrix with non-negative 
entries and column sums equal to one. Each column of such a matrix is 
a probability vector with a corresponding Shannon entropy. The minimal 
output entropy is easily seen to be the smallest of these entropies, in other 
words, the minimizer is the degenerate measure corresponding to the mini- 
mizing column. Another well-known case is that of Bosonic Gaussian chan- 
nels. It was conjectured for a long time that the minimizer is a Gaussian 
state. This conjecture was recently proven in [19]. Also for Davies maps 
there are natural candidates for a minimizer: the eigenstates of the system 
Hamiltonian. This turns out to be false even in the qubit case, at least for 
maps that are sufficiently close to the identity. After a sufficiently long time 
the Davies map will generically be close to the projection on the equilibrium 
state and therefore behave more or less like a classical stochastic map. One 
indeed observes that in such a situation the minimizer is an eigenstate of 
the Hamiltonian. 

Let us take the input state 



V ) , < n < 1, and < v < a/M 1 - ju)> 
1 /iy 



Davies map 

where b = ap/(l — p) like in f)3. 12|) and (|3. 13j) . one gets 
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— fi) . After actin 
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— a 


b ^ 









c 









c 
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1-bJ 





p 



(1 — a)fi + 6(1 — /i) cu 

cv a/i + (1 — 6)(1 — /i) 



The eigenvalues £± of pi are 

£± = ^ ( 1± V 1 - 2fo + 2 (- 1 + a + % + 4c2l/2 )- 

The entropy reaches its minimum when the difference between the two eigen- 
values is maximal. Therefore we need to find the maximum of the expression 
under square root substituting z/ 2 by /x(l — /x). The second derivative of this 
expression is equal to 8(1 — a — b) 2 — c 2 . A number cases have to be distin- 
guished: 
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1. If c 2 < (1 - a - b) 2 then the minimizer is fj. = which is the ground 
state of the system Hamiltonian and the corresponding output state 
is diagonal 

'1-6 N 



*(/W)=[ Q b j. (3.19) 

2. If (1 — a — b) 2 < c 2 < (1 — a — b) (1 — 26) the minimizer is again 
[i = 0. In this case the function attains its maximum outside the set 
of admissible values of \i and therefore fi is an endpoint of (0, 1). 

3. If (1 - a - 6)(1 - 26) < c 2 the the minimizer is determined by 

_ (a + 6- 1)(26- 1) -c 2 
M0 ~ 2(a + 6 - l) 2 - 2c 2 ' 

In this case it is a true superposition of ground and exited states of 
the system Hamiltonian. The eigenvalues of the output state are 




c 2 (c 2 + o(2 - 46) + 26- 1) 
c 2 - (-a -6 + 1 



e± = o(i±a/ \. 2 ; : z i ^2 )• ( 3 - 2 °) 



The entropy is easily computed from the formulas (|3.19p and (13.2QH . The 
condition c 2 < (1 — a — 6) 2 written in terms of decay constants translates to 

n < r 3 . 



4 Davies maps on a three-level system 

As discussed at the end of Section [231 a generic Davies map $ stochastically 
mixes the diagonal components of a density matrix and scales down the 
off-diagonal entries. The full Davies superoperator acting on a three-level 
quantum system is therefore represented by a matrix 





— F21 — F31 
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-F32 











1 - ^13 - 

(4.21) 



^23/ 
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The stochastic matrix 



(1 — F21 — F31 F\2 i*i3 

F 21 l-F 12 -F 32 F 23 j (1.22) 

F31 -F32 1 — F\s — F 2 3 y 

satisfies the classical detailed balance condition with respect to the proba- 
bility vector (pi,p 2 ,P3) where the pi are the eigenvalues of pp, determined 
by the energy levels of the system Hamiltonian and by the temperature 

/•:,/;, /•),/>;. i,j = 1,2,3. (4.23) 

As in the qubit case we assume that the eigenvalues of the system Hamil- 
tonian are strictly ordered: e± > e 2 > e%, put differently p\ < p 2 < P3- 
Assuming that all entries of F are strictly positive, the conditions (14. 23ft are 
satisfiable if and only if 

F\ 2 < F 2 i, F13 < F31, F 2 3 < F32, and F± 2 F 2 3 F31 = F13 F32 F21, 

one of the inequalities being redundant. The stationary vector is given by 



Pi 



Fl 2 F\3 + F12F31 + ^13^21 



and the other components p 2 and P3 are obtained by cyclically permuting 
this expression. 

Another way to look at the stochastic map F is to say that it is a special 
case of a quantum map: it mixes the diagonal elements of a density matrix 
according to F and destroys the off-diagonals. Complete positivity is auto- 
matically satisfied because the Choi matrix of this map is diagonal with the 
entries of F on the diagonal. 

Furthermore, if $ is a Davies map then it lies on a Davies semi-group and this 
implies that F lies on a classical detailed balance semi-group of stochastic 
maps 

F = exp(G) 

with G a classical detailed balance generator and t = 1. In particular this 
means that 

G12 

G = I G21 — G\ 2 — G32 

G32 
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with 

Gij Pj = GjiPi, i,j = 1,2,3. 

As argued at the end of Section 12.31 detailed balance generators have real 
eigenvalues. A property which we have to check in order to ensure the 
existence of a semi-group F = exp(G) is positivity of the eigenvalues of 
F because only in this case there exist a real logarithm of F. The three 
eigenvalues of F in decreasing order are 

a(F) = {1,A + B,A-B}, (4.24) 

where 

A = \ (Tr F - 1) and B = \ y / 2Tri ?2 — (Tr F) 2 + 2 Tr F - 3. 
The existence of a logarithm therefore requires the positivity conditions 

< A and B < A, 

these are equivalent with 

l<TrFand Tr F 2 + 2Tr F < 2 + (TrF) 2 . 



The last criterion which has to be satisfied is that G = log(-F) is a generator 
which means that its off-diagonal elements must be positive. It is possible 
to compute log(F) straightforwardly, see the appendix. The last condition 
that guarantees that F belongs to a semi-group reads 

G12 > 0, G 23 > 0, and G 31 > 0. (4.25) 

These conditions are rather complicated. There are three of them relating 
the spectrum (|4.24|) of F with its entries. We just write down a single one 
explicitly derived in the appendix: 

Vl {l - A + B)\og{A + B) < y 2 (l - A - B) \og{A - B) (4.26) 

where 

Vl := 2B - F 12 - F 21 + F 13 - F 31 + F 23 - F 32 + 2 %^1 and 
V2 '■= 4 — 2B — F\ 2 — F 2 \ + F\ 3 — F 3 \ + F 23 — F 32 + 2 — . 

-^21 
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The other two conditions are obtained by cyclically permuting that of above. 

We now return to the full quantum map <£. Clearly the Xj in (|4.21j) must 
be positive as they have to be exponentials of real numbers. We write 
Tj = — logAj. An additional condition should be imposed in the quantum 
case. The generator of <3? has to be conditionally complete positive which 
means that 



Together with (|4.25p the above relations guarantee that $ is a Davies map. 
We assumed that the Bohr spectrum of the system is non-degenerate, that 
is that there are no degeneracies in the non-zero differences of energy levels 
of the Hamiltonian. If it is not the case then blocks appear on the diagonal 
of 

5 Conclusion 

We studied a class of quantum operations called Davies maps. They satisfy 
the quantum detailed balance condition and the semi-group property = 
exp(Qt). Moreover, the unitary and the dissipative parts of the generator 
Q commute. This property implies that the dynamics of populations, the 
diagonal elements of a density matrix p, separates from the evolution of 
coherences, the off-diagonal elements of p. 

For a single qubit the complete characterization of the set of Davies maps 
has been worked out. This three dimensional set of maps is determined by 
parameters a and c subject to the conditions (|3.13p at a given temperature 
T. These parameters determine the lengths of the axes of the ellipsoid 




This translates into the conditions 



n + r 2 + r 3 > G12 + G 23 + G 3 i + Gn + G 32 + G w 



and 



t\ + r\ + rf + 2(n(G 21 + G 31 ) + r 2 (G 32 + G 12 ) + r 3 (G 3i + G 12 )) 
< 2(rir 2 + r 2 r 3 + r 3 n) + ((G 2 i + G 3 i)(G 32 + Gi 2 ) 
+ (G 32 + G 12 )(G 13 + G 23 ) + (Gi 3 + G 23 )(G 21 + G 31 )) 
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obtained by applying the Davies map $ to the Bloch sphere. Thus they 
are related to both relaxation rates t\ = T2 and T3 often used to describe a 
non-unitary one-qubit map in quantum optics. The constraints (|3.17p which 
determine the set of Davies maps are found to be equivalent with conditions 
for a master equation of the Lindblad type. 

Alternatively, one may parametrize the set of one-qubit Davies maps by the 
following three numbers: the vertical shift K3 of the centre of the Bloch ball, 
see Fig. HI the volume of the image of the ball and its asymmetry parameter. 
It is easily seen that the first parameter depends of the temperature, while 
the two others are functions of a and c. 

So far we restricted our attention to the low dimensional cases A = 2 and 
N = 3 but some remarks can be made on general A- level Davies maps. 
Consider such a system described by a Hamiltonian with a non-degenerate 
spectrum. A Davies channel that satisfies all conditions specified in Sec. 12.31 
is described by a block matrix of size A 2 . The outer block is a stochastic 
matrix F which describes the transfer of populations. This matrix obeys 
the classical detailed balance condition, see (|4.2ip and (|4.22p written for 
iV = 3. The number of parameters needed to describe F is (A 2 + A)/2 — 1: 
the number of elements on the diagonal of F and above it minus one due 
to the trace-preserving condition. There are also (iV 2 — A)/2 elements of 
$ which describe the evolution of quantum coherences, the parameters Aj 
in (I4.2ip . Together this gives d = A 2 — 1 parameters necessary to specify 
a Davies map acting on an A-level system. This is exactly the number of 
parameters required to characterize a density matrix of this size. Thus the 
Davies maps forms a measure zero subset of all stochastic maps, a set of 
dimension D = A 4 - A 2 = N 2 d. 

We computed the minimal output entropy of a one-qubit Davies map and 
found that in the quantum regime, close to the identity map, the minimizer is 
not equal to the obvious candidate: the eigenstate of the system Hamiltonian 
which remains as pure as possible under the stochastic mixing of diagonal 
elements. It is then interesting to ask for the additivity of minimal output 
entropy for this class of maps. Independently whether this property holds 
true, we believe that Davies maps form a physically relevant class of quantum 
operations, an interesting subject of further study. 
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A Derivation of inequality (14.261) 



Let X be a 3 x 3 complex matrix and assume for simplicity that it has non- 
degenerate eigenvalues £i, £21 and £3. This implies that X is diagonalizable. 
Let g be a complex function defined on a domain containing the eigenvalues 
of X, then 

g(X) = got + giX + g 2 X 2 

where go, g\, and g 2 are complex functions that only depend on the eigen- 
values of X. More explicitly 

go = -77 — — ft #(6) + c y cl - P erm - 

(€1 - ?2M3 - SiJ 

91 = u FYTc F^3{Ci) + cycl. perm. 

(£1 - £2 J (£3 - £ij 

92 = -77 -r-rrz T\9(£,l) + cycl. perm. 

(£1 - ?2j(?3 - £1) 

We can use this to compute the logarithm of the stochastic matrix F in (I4.22|) . 
Using the parametrization (|4.24p for the eigenvalues of F we obtain expres- 
sion ([426]). 
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